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ABSTRACT 

We study the spectral and temporal behavior of X-ray flares from the active M-dwarf EV Lac in 200 
ks of exposure with the Chandra /YWTGS. We derive flare parameters by fitting an empirical function 
which characterizes the amplitude, shape, and scale. The flares range from very short (< 1 ks) to 
long (~ IC* s) duration events with a range of shapes and amplitudes for all durations. We extract 
spectra for composite flares to study their mean evolution and to compare flares of different lengths. 
Evolution of spectral features in the density-temperature plane shows probable sustained heating. 
The short flares are significantly hotter than the longer flares. We determined an upper limit to the 
Fe K fluorescent flux, the best fit value being close to what is expected for compact loops. 

Subject headings: X-rays; stars; spectra; stars:individual EV Lac 
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1. INTRODUCTION 

Coronal activity is ubiquitous among late type stars of 
all classes: pre-main sequence, main sequence, evolved 
binaries, and even some single giants. Activity on main 
sequence stars is generally correlated with rotation rate 
except at the shortest periods where saturation occurs. 
Flaring is common among the coronally active stars 
(as defined by high X-ray band luminosity relative to 
the bolometric luminosity, log(Lx/-^boi) ^ —4 to —3). 
Nearly every sufficiently long (many kiloseconds) X-ray 
observation of coronal sources shows flares, as defined 
by a rapid increase in flux accompanied by a harden- 
ing of the spectrum and subsequent decay and softening. 
Flares are the most dynamic aspects of coronal activity 
and are possibly a significant source of coronal heating. 
From Solar studies, we know that flare mechanisms in- 
volve magnetic field reconnection, particle beams, chro- 
mospheric evaporation, rapid bulk flows, mass ejection, 
and heating of plasma confined in loops. 

Many stellar activity studies have attempted to avoid 
flares in order to determine quiescent coronal properties 
such as emission measure distributions, elemental abun- 
dances, loop heights, and geometric distribution of ac- 
tive regions. Since flares are a phenomenon of magnetic 
reconnection thought to occur on small spatial scales, 
modeling the dynamic behavior allows us to constrain 
loop properties in ways that cannot be done from analy- 
sis of quiescent coronae that necessarily require a spatial 
and temporal average over a large ensemble of coronal 
structures. Use of rotational modulation of flux and ve- 
locity has become a common and important technique 
(Doppler imaging, and related methods) for mapping 
the stable and non-uniform structures of stellar activ- 
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ity. Flare modeling has the potential to become as 
important for determining the properties of transient 
loop structures and their ene rgetics (e.g. iReald 120071 : 
lAschwanden fc Tsiklaurl[2009l ) . 

Here we exploit the flaring behavior in EV Lac — one of 
the brightest, most reliable flaring sources available — to 
obtain spectra and temporal profiles in flares from Chan- 
dra X-ray Observatory (CXO) High Resolution Trans- 
mission Grating Spectrometer (HETGS) observations. 
The 200 ks exposure, flare frequency, and amplitudes are 
sufficient to provide spectral diagnostics in flares if we 
combine multiple events to examine the mean properties 
of emission line and continuum evolution. In addition, 
we study the distribution of flare temporal morphologies. 

Even though EV Lac has a high flare rate, we are still 
forced to work with some mean quantities by combining 
data from multiple, possibly physically different flares. 
Nevertheless, the results are interesting and are neces- 
sary groundwork for guiding future studies. In this pa- 
per, we characterize the flares and present some simple 
diagnostics from different flare states. We also model 
Fe K fluorescence. Detailed hydrodynamic models will 
be applied in future work. 

1.1. Characteristics of EV Lac 

EV Lac is a nearby (5 pc) dM3.5e (O.SSAf©, O.36i?0) 
flare star with a photometric period of 4.4 days. It 
is among the X-ray brigh test of single dMe flare stars 
(jRobrade fc Schmittl[2005[) . having a mean Lx ^ 3 — 5 x 
10^® ergs cm~^ s~^. It is of particular interest for this 
study for these specific reasons: 

• EV Lac has consistently shown frequent and 
strong flares whenever observed in X-ra ys, as has been 
demo nstrated wi th ROSAT a i id SA X (jSciortino et al.l 
[1991 ■ ASCA (iFavata et all [2000h. X MM I Newton 
(IMitra-Kra ev et all 120051: IRobrade fc SchmittI 
I2006D. Chandra (lOsten et all I2Q05D . and Suzaku 
( Laming fc Hwan^l20Q9D . Models imply that the flaring 
structures are of a compact nature. Such a compact 
geometry is more efficient for producing Fe K fluores- 
cent emissions due to the larger solid angle seen by 
the photosphere, and also because lower loop heights 
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have a greater the yield of Fc K photons from the 
fluorescing photosphere — more hard X-ray photons 
enter at small an gle and have a lower optical depth for 
escape (lBail[T979l ). The X -ray flare frequency is about 
0.4 hr~^ (see Figure [1]), similar to previously quoted 
rates of 0.2 — 0.3 hr~^ for X-ray, UV, and U-b and flares 
(jOsten et all 120051 iMitra-Kraev et alj I2005D . Some 
of the flare s are very short (< 1 ks; Figure [T] and 
lOsten et al.l (I2005D ). The short duration attests to a 
compact size since larger volumes have longer cooling 
times due to their lower density which reduces the 
radiative and conductive loss rates. With the addition of 
the new HETG observations, we have identified 25 flares 
in 200 ks (Figure [T|) , with a broad range in amplitude 
and timescales. 

• EV Lac's O VII triplet ratio showed a dens i ty of about 
6 X 10^° cm-3 (jTesta. Drake fc Peres! 120041 iNess et al.l 
[200l . The ability to measure density provides strong 
constraints on the emitting volume. 

• EV Lac is one of 3 in a sample of 22 coronally ac- 
tive stars which showed presence of opacity as deter- 
mined &ran_thfiratiosof H-like Lya to Ly/3 resonance 
lines (jTesta et al.ll2007al ). Measurement of opacity can 
also provide constr aints on emitting structure geometry. 
iTesta et "all ()2007a|) derived a compact coronal height of 
about 0.06 stellar radii. 

• EV Lac has a strong and variable magnetic 
field (3kG) with a bout a 50% surface filhng f actor 
(jJohns-Kruh fc Valc nti 199a iPhan-Bao et al.ll2006D. and 
it w as both poloidal and asymmetric (jPonati et al.l 
I2008D . Strong fields may be necessary for forming and 
maintaining compact loops which give rise to the short, 
energetic hares. The star is near the expected mass 
boundary for transition from an a — w dynamo to a tur- 
bulent dynamo and so could be an important case theo- 
retically. 

2. Chandra/UETGS OBSERVATIONS AND DATA 
PROCESSING 

Chandra has observed EV Lac twice with the HETGS, 
once in 2001 (Observation Identifier 1885) for 100 ks, 
and once in 2009 (ObsID 10679, as part of the 
HETG Guaranteed Time Observation program) for 
97 ks. T he spectrometer and obse rvatory are de- 
scribe d by iCanizares et al.l (|2005| ) and IWeisskopf et al.l 
(|2002( ). respectively. The observational data were 
processed with the standard software suite, CIAO 
(fFruscion c et al.l 120061 ) to filter, transform, bin spec- 
tra, and construct the observation-dependent responses. 
Thes e CIAO programs were in turn driven by t he TG- 
Cat (jMitschang. Huenemoerder fc Nichols l2010f ) repro- 
cessing script^ which automates the CIAO processes 
into an end-to-end pipeline. This was especially useful for 
extraction of hundreds of spectra and responses in time- 
filtered intervals (see Q ■ Light curves were binned from 
the event files using the "ACIS Gratings Light Curve" 
package (aglc0), which bins counts and count rate light 
curves over multiple detector chips, grating orders, grat- 
ing types, and wavelength regions of the dispersed grat- 
ing spectrum. 
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3. FLARE LIGHT CURVE FITTING & STATISTICS 

The WcibuU distribution is a convenient function for 
empirically fitting the shapes of the fiares, whether im- 
pulsive or gradual. The normalized distribution is 

/(p;a,s)=(f)p(''-i)e-P" (1) 
{t-to)/s (2) 

in which a is a shape parameter (a > 0), the decay scale 
(or width of the distribution) is specified by s (s > 0), 
the offset is given by to (to > 0); the independent time 
coordinate is t (t > to). Shapes a < 1 arc exponential- 
like — there is no resolved rise. A shape a > 1 is more 
Gaussian-like, becoming more symmetric as a increases. 
For a < 1, the function falls faster than an exponential, 
and for a > 1, it falls somewhat slower than an expo- 
nential, up to a factor of about 1.5, as defined by the 
e-folding time from the peak rate. 

Light curves were binned to 100 s using events from 
positive and negative first orders, HEG and MEG, from 
1.5-25 A. These uniformly-binned light curves were fit 
iteratively using a sum of Weibull distributions scaled by 
an amplitude parameter, and including a constant basal 
count rate; in other words, we fit Af(p; a, s) + Ro in 
which A is the flare area in counts, Ro the baseline rate 
in ctss~^ for the observation, and / is given by equa- 
tions [lEl We flrst dcflned model components for obvi- 
ous flares, fltting the region of significantly overlapping 
flares, and then added components as required to flat- 
ten the residuals. Since we physically expect single-loop 
flares undergoing radiative decay to be exponential in 
time, or likely prolonged by sustained heating, we con- 
strained the shape, a, to be > 1 so that the decay is 
exponential or slower. 

Table [T] gives the fiarc parameters. The first column, 
n, is an arbitrary index used as a unique identifier for 
each flare (with no signiflcancc to the order). The sec- 
ond column, area, is the overall normalization, or ampli- 
tude, of the flare in counts, integrated over 1-25 A. The 
next three columns (3-5) give the Weibull distrubution's 
parameters, to, the flare start time (measured from the 
start of the observation), the shape parameter, a, and 
the scale parameter, s. The values in parentheses are 
the 90% confidence intervals of the parameters. 

The last three columns (6-8) are not independent pa- 
rameters, but are given as an alternate and more familiar 
characterization of the flares. They are the e-folding time 
from the peak, ti , the lag from the start to the peak, T2 , 
and the count rate at the peak, r. The horizontal line 
midway down the table separates the 2001 observation 
(n > 100) from that of 2009. Figure [2] shows the shapes, 
scales, and areas graphically. The baseline count rates, 
Ro were slightly different for the two observations, being 
0.151 (0.004) in 2001, and 0.092 (0.002) in 2009 (uncer- 
tainties are for 90% confidence). 

4. FLARE SPECTRAL PROFILE EXTRACTION & FITTING 

In order to get enough counts in different flare states, 
we need to extract spectra as a function of time during 
flares, then form composite spectra from states deemed 
similar. Time ranges and number of spectra to extract 
were identified manually from the light curves, using 
higher time resolution during intervals of more rapid 
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Figure 1. These arc the light curves from the two HETGS observations, for HEG and MEG 1 — 25 A events. The top shows the 2001 
observation (ObsID 1885), and bottom 2009 (ObsID 10679). It is clear that flares are frequent; 25 were identified and fit with an empirical 
function, for which we show the model (overlayed smooth, light-colored curve) and residuals (lower panels). 



Table 1 

Flare Light Curve Fit Parameters 



n 


area 
[cts] 




to 








a 




s 

N 




r2" 


[cts/s] 


(1) 


(2) 




(3) 








(4) 




(5) 


(6) 


(7) 


(8) 


101 


595 (494, 713) 




(0, 524) 


1 


89 


(1.46, 


2.47) 


6474 (5691, 7857) 


9818 


4346 


0.08 


111 


234 (183, 292) 


15685 


(15685, 


15785) 


1 


40 


(1.04, 


1.77) 


1383 (1047, 1833) 


2116 


563 


0.12 


112 


226 (170, 281) 


19101 


(19101, 


19202) 


1 


32 


(1.00, 


1.75) 


1715 (1304, 2144) 


2586 


588 


0.10 


102 


1533 (1452, 1614) 


40861 


(40747, 


40904) 


1 


85 


(1.72, 


2.02) 


2181 (2068, 2316) 


3320 


1432 


0.57 


103 


854 (802, 905) 


45726 


(45721, 


45726) 


2 


13 


(2.00, 


2.30) 


239 (230, 250) 


354 


177 


3.18 


104 


184 (138, 232) 


49645 


(49486, 


49745) 


1 


47 


(1.02, 


2.81) 


698 (508, 866) 


1076 


321 


0.20 


105 


1154 (1073, 1237) 


51548 


(51548, 


51926) 


1 


89 


(1.55, 


2.05) 


3752 (3358, 3954) 


5689 


2517 


0.25 


106 


675 (554, 802) 


60934 


(60898, 


60998) 


1 


00 


(1.00, 


1.13) 


3100 (2322, 4115) 


3104 





0.22 


107 


666 (586, 749) 


66725 


(66642, 


66825) 


1 


53 


(1.33, 


1.72) 


1659 (1462, 1840) 


2563 


834 


0.30 


108 


1753 (1669, 1840) 


76370 


(76364, 


76370) 


1 


45 


(1.38, 


1.52) 


819 (777, 868) 


1260 


366 


1.58 


109 


4639 (4492, 4786) 


78078 


(78065, 


78078) 


1 


19 


(1.15, 


1.22) 


3264 (3140, 3397) 


4627 


682 


1.08 


110 


3089 (2928, 3260) 


87418 


(87255, 


87422) 


1 


38 


(1.31, 


1.47) 


7105 (6681, 7611) 


10839 


2776 


0.32 


9 


7742 (7436, 9336) 




1305 (0, 1371) 


1 


00 


(1.00, 


1.01) 


8071 (7652, 8529) 


8084 





0.96 


8 


1052 (725. 1292) 


14639 


(14538, 


14639) 


1 


01 


(1.00, 


1.36) 


2738 (1505, 3429) 


2912 


33 


0.36 


10 


200 (123, 450) 


17654 


(17300, 


18358) 


2 


45 


(1.00, 


4.81) 


1895 (1082, 2855) 


2724 


1528 


0.10 


4 


537 (454, 645) 


34448 


(34261, 


34540) 


1 


19 


(1.01, 


1.44) 


3501 (2853, 4718) 


4988 


761 


0.12 


1 


1291 (1140, 1416) 


42782 


(42664, 


42783) 


1 


24 


(1.15, 


1.44) 


3804 (3366, 4469) 


5575 


1024 


0.25 


12 


1154 (1035, 1296) 


42969 


(42891, 


42984) 


1 


71 


(1.56, 


1.99) 


890 (827, 994) 


1369 


533 


1.02 


2 


760 (702, 818) 


51628 


(51601, 


51628) 


1 


24 


(1.16, 


1.33) 


616 (563, 674) 


904 


167 


0.92 


3 


166 (127, 232) 


55447 


(55347, 


55548) 


1 


16 


(1.00, 


1.51) 


1202 (898, 2009) 


1664 


212 


0.11 


13 


547 (428, 691) 


60975 


(60900, 


61478) 


1 


42 


(1.09, 


1.71) 


6547 (5303, 8475) 


10046 


2787 


0.06 


7 


71 (50, 181) 


64091 


(63982, 


64091) 


1 


27 


(1.00, 


3.33) 


142 (100, 225) 


210 


42 


0.37 


6 


406 (304, 511) 


68413 


(67939, 


68814) 


1 


81 


(1.46, 


2.47) 


3918 (3294, 4760) 


5981 


2512 


0.08 


5 


436 (397, 475) 


79369 


(79354, 


79369) 


1 


57 


(1.45, 


1.72) 


424 (390, 460) 


655 


224 


0.78 


11 


2246 (1012, 63190) 


88716 


(88315, 


88918) 


1 


12 


(1.00, 


1.28) 


25389 (10008, 278646) 


33882 


3459 


0.07 



Note. — The central horizontal line separates the 2001 and 2009 observations. Parameters to, a, and s are as defined by 
Equations lH2l and the area is the flare amplitude. The index n is an arbitrary idcntiflcr (with no significance to the order); it 
is used to mark the fiares displayed in Figure[3] The vertical line separates the fitted parameters from some derived parameters 
which may be more intuitive. Values in parentheses are 90% confidence limits. 

^ ri (the e-folding time), T2 (time from to to the peak rate), r (the peak rate), are not unique parameters, but arc derived from 
the preceding WeibuU distribution's parameters. 



change. Given a list of start times, stop times, and num- 
ber of spectra, good-time- interval (GTI) tables were con- 
structed which were then used to filter the events and 
extract spectra. We extracted 121 spectra for the 2001 
observation, and 118 for 2009. Processing was automated 
using the tgcat scripts so that each spectrum also has 
the appropriate response for its time interval. 

Given the extractions, the flare phase of each time bin 
was then defined as "low" , "rise" , "peak" , "decay high" , 
"decay low", or "ignore". Of the 249 total spectra ex- 
tracted (each comprised of HEG and MEG positive and 



negative first orders), 109 were assigned to the noticed 
groups, the ignored remainder being of ambiguous phase 
due to overlapping flares of comparable magnitude. Ta- 
bic [2] gives some summary statistics on the selections. 
Since the rises are often very sharp, we have the shortest 
exposure and fewest counts in this group. The longest 
time bin, for the low phase, was 17 ks, while during flares 
we used bins as short as 300 s. In Figure|3]we show detail 
of some of the flare phase selections. 

The spectral analysis fitting was done by loading 
the individual spectral counts histograms and their 
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Shape 

Figure 2. Flare fit parameters. These specify tlie model used in 
Figure [T] and listed in Table [T] A shape (a) of 1 means exponential 
decay, while larger shapes have resolved rises. The scale (s) is 
approximately the e-folding time from the peak rate. Circles' radii 
are proportional to the log of the flare counts, grouped in quartiles 
(the inset key gives actual ranges in counts). There is no obvious 
correlation among these parameters. 
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the index, n, as given in Table [Tl 



associated respons e files into the ISIS (IHouckl [200l 
iHouck fc Denicolal 120001 ) analysis system where they 
could then be associated and combined dynamically dur- 
ing fitting. Spectral grouping was also done dynamically 
as appropriate to obtain sufiicicnt statistics per wave- 
length bin. 

5. FLARE MODELING 

We study the flares by applying different model 
approaches to derive diagnostics at different levels, 
using Solar flare models as a basis. The numer- 
ical models have become very sophisticated; codes 



Table 2 

Flare Phase Information 



Phase 


N 


^cxp 

[ks] 


Counts 


low 


13 


50 


7235 


rise 


15 


9 


3036 


peak 


21 


10 


8133 


decay-high 


26 


18 


10585 


decay-low 


34 


22 


7369 



A'^ is the number of spectra extracted in 
each flare phase. 




4x10-' 



5x10" 6x10* 
Time [s] 



Figure 4. Here we show detailed evolution for flares 104-106, 
events with resolved rise and decay, no resolved rise, or only the 
unresolved peak. We define the flare "phases" as the rise (trian- 
gles), peaks (stars), decay-high (squares), decay- low (circles), and 
low (crosses). 

and methods have been tested in detail against 
spati ally resolved solar observations (see for exam- 
ples. lRealell2Q07tlKlimchuk. Patsourakos fc Cargilil [200l 
lAschwanden fc Tsiklaurill2009f ). In application of these 
models to EV Lac flares coupled with spectral diagnos- 
tics, we can potentially constrain flare loop conditions in 
this M-dwarf and relate it to the Sun and other stars. 

Many solar X-ray flares are characterized by simple 
rise plus decay light curves and typically involve local- 
ized loop structures, where the plasma is confined by 
the coronal magnetic field. It is then possible to de- 
scribe the flaring plasma as a compressible fluid which 
moves and transports energy along the magnetic field 
lines and use time-dependent hydrodynamic models. In 
these conditions, it has been shown that the flare decay 
time scales with the length of the flaring loop (|Serio et al.l 
[T99l[ ) and the presence of significant heati ng in the decay 
can make the decay longer than expected (jJakimiec et al.l 
[199 2: Svl wester et al. 1993). Such heating can be diag- 
nosed from the analysis of the fiare path in the density- 
temperature diagram: the flare decay path becomes shal- 
lower. From numerical modeling, scaling laws have been 
derived to estimate t he loop length af ter correcting for 
the effect of heating (jReale et al.lll997D : Lg = CrxV^V, 
where Lg is the loop half length (10^ cm), the de- 
cay time of the X-ray light curve in seconds, Ty is the 
maximum flare temperature (10^ K), and C a propor- 
tionality factor. Such a diagnostic tool is well estab- 
lished and has been commonly applied to analyze stel- 
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kr flares (iReale fc Micel alll998l: [ Fayata fc Scliniittlll999l: 
Maggio et alJ 120001: "iStelzer et all l20Q2t iBriggs fc Pvd 
2Q03t iPillitteri et al.|[20b5i : iFavata et al.ll20Q5D. Withou t 
heating, the scale factor C = 1/120 (|Serio et al1ll991f ). 
but the scaling law can greatly overestimate the loop 
length. More recently, new diagnostic tools have been 
developed to obtain information on the flaring pla. sma 
and structures from the flare rise phase (|Realell2007[) . It 
has been also shown that flares with more complex light 
curves, e. g:. with m ultiple peaks, can involve more coro- 
nal loops (jReale et al. 2004). 

We are limited by statistics to composites of hetero- 
geneous flare types. Nevertheless, the behavior agrees 
in general with theoretical expectations. Since no sin- 
gle flare is strong enough for this analysis, we combined 
spectra from several flares and provide diagnostics for 
both lines and continua. One theoretically interesting 
diagnostic is the evolution of temperature and density 
(or its proxy, V EM). The helium-like to hydrogen-like 
line ratios are a strong function of temperature. From 
a portion of the spectrum which includes a line pair, 
wc can derive a characteristic temperature and emis- 
sion measure from a relatively simple, single-temperature 
APED model fit. Figure [S] shows the Si Xiv-Si xiii re- 
gion. We have derived temperatures from fits to sev- 




eral H-like to He-like line pairs and a 3 — 6 A continuum 
band using APED emissivities. Figure [6] shows the re- 
sults for many flares in the existing data, grouped by 
flare phase as given in Table [2] These agree q ualita- 
tively with paths expected from single-loop theory (jReald 
I2007D . The initial low-phase is the left- most point for 
each feature, and flare phase (or time, if it were a sin- 
gle flare event) proceeds (generall y) clockwise around the 
loop. This is similar to results of lTesta "eFall (|2007bD on 
a single, large flare in the active giant, HR 9024. Fig- 
ure [5] also shows two arbitrarily placed lines, one of slope 
2, and the other of slope 1/2. The former is charac- 
teristic of radiative and conductive decay, and the lat- 
ter of quasi-steady-state conditions in which the heating 
changes slowly enough that the temperature and den- 
sity can obtain th eir equilibrium values in a hydro static 
magnetic "RTV" (|Rosner. Tucker fc Vaianal[T97l loop. 
( Fo r details on su c h hyd r odyna,niic rn o del trajectories, 
seelJakimiec et al.l (|1992[ ): iReali ()2007[ ): ISvlwester et al.l 
(|1993[) .) The measurements shown for EV Lac generally 
have slopes between these limits, implying that our flare 
composite has some sustained heating. The fact that 
some of the lines cross (e.g., for Si and Mg) likely means 
that we have mixed flares of very different physical char- 
acteristics, and that we require flner grouping by flare 
types (and consequently, better statistics). 

The temperatures at identical phases differ between 
spectral features because we are sampling from plasma 
with a distribution of temperatures, and the features 
have different emissivity dependence upon temperature. 
In Figure [51 if we were to connect the points at identical 
phases, then we would have a crude emission measure 
distribution with temperature (in the square root, and 
vertically oriented) for each flare phase. 

The naive scaling law for the range in peak tempera- 
tures seen (Ty ~ 0.6—3.0) and broad range in decay scales 
(200 — 30000 s) would imply a huge range in ig, of two 
orders of magnitude, ~ 1-100. However, we know that 
we have a heterogeneous mix and that there is sustained 
heating, so the upper value is clearly a gross overestimate 
(for reference, the stellar radius in these units is 25). If 
we consider only the shortest flares — which are more 
likely to occur in one or few loops — then the range is a 
more plausible Lg ~ 0.1 — 1. 

Further progress will likely require us to obtain bet- 
ter statistics in groups of similar flares, in order to con- 
strain loop parameters with hydrodynamic models. For 
instance, due to limited time resolution, we are not 
certain that the temperature and density peak at the 
same flare phase, as they appear to do in Figure [51 If 
we had better temporally-resolved spectra for several 
similar flares (shape and duration), then we could ap- 
ply detailed hydrodynamic modeling by solving accu- 
rately the time-dependent hydrodynamic plasma equa- 
tions with the aid of the Palermo-Har vard numerical code 
(jPeres et al.lll982l: iBetta et al.lll997[) . This approach al- 
lows one to obtain deeper insight into the flaring plasma 
and of the heating details, such as the detailed ther- 
mal structure and its evolution, an d information about 
Reale et al.l [200l iFavata et all 120051 : 



the loop aspect 
iTesta et al.l 120071^ 



Figure 5. This figure shows the evolution of the He-hkc Si xiv 
(6.18 A) to Si xiu (6.65 A) fine fluxes through a composite flare. 
Note that the y-axes have very different scales. 



Such will be attempted in future 
work, and hopefully with a larger dataset which supports 
grouping of more like events. 
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Figure 6. Evolution of temperature and density proxy using sev- 
eral line pairs or a continuum band for the composite flare phase 
selections. The line flux and temperature diagram traces the evo- 
lution of a flare, using composite spectra from different phases of 
many flares. The 5 points on each curve correspond to the dif- 
ferent flare phases, as defined in Table (2] starting from the low 
phase at the left-most point in each case and proceeding through 
"rise" , "peak" , "decay-high" , and "decay-low" . When combined 
with time-dependent hydrodynamical flare models, the slopes in 
the different phases can be related to flare loop sizes and heating 
functions. We show two limits as gray lines — a slope of 2 is char- 
acteristic of radiative and conductive decay, and a slope of 1/2 of 
quasi-steady-state decay. Errorbars are the 90% confidence limits 
from the single-temporaturc APED fits to each spectral region. 

5.1. Spectral Comparison of Different States 

Wc can empirically compare spectra for the short and 
long flares to each other and to the low, or "quiescent" 
phase. We have extracted spectra integrated over the 
flare duration for short flares (flare numbers 2, 5, 7, 12, 
103, 104 and 108; see Figure [3]), long flares (numbers 
4, 6, 8, 9, 10, 102, 105, 109, 110), and for low phases. 
In these groups, we respectively have exposures of 7ks, 
42 ks, and 50 ks. To characterize the spectra, we have fit 
a 3-component APED model to each. Such a model is not 
as definitive as a line-based emission measure and abun- 
dance reconstruction. The temperature represents some 
average value since the line features can change dramati- 
cally over temperature differences of order of a factor of 2, 
which a 3-temperature model cannot fairly represent over 
the broad range of plasma temperatures. Also, the abun- 
dances are largely another line-strength parameter which 
can compensate for off-nominal temperatures as well as 
real abundance values; we do not consider abundance to 
be an interesting parameter for this purpose (determina- 
tion of meaningful abundances requires line-based emis- 
sion measure reconstruction). Nevertheless, the model 
is sufficient to show the primary differences between the 
phases. The short wavelength continuum (< 6 A) con- 
strains the highest temperature component, while the 
lines, primarily the H-like and He-like series, constrain 
the lower temperature components. Some lines, such as 
Fe XXIV can also have a significant effect at the high 
temperatures of flares. 

In Figure [71 we show the three different spectra for 
different wavelength intervals. The upper (at the shorter 
wavelengths) curve (blue) shows the short flares' peak 
phase. It is somewhat noisier than the other two curves 



due to its lower exposure, but is systematically stronger 
and has features characteristic of high temperature which 
are weaker or not present in the long flare (red middle 
curve) or quiescent (black lowest curve) states. Fe xxv 
and the S xv-S xvi series are prominent during the short 
flares (upper panel), as are other high ionization states 
of iron, such as Fe xxiii and Fe xxiv (second panel). At 
longer wavelengths (lower two panels), where low tem- 
perature features dominate, changes are less pronounced. 

In Figure [8] we show the temperature versus emission 
measure from the three-component APED fits. It is clear 
that the shortest flares have much greater weight in the 
hottest plasma than the long flares. The long flares do 
have some high temperature plasma, but it does not 
dominate the emission measure. Their greater weight in 
cooler plasma during the decay implies that the emission 
originates from a multi-loop fl aring system with simul- 
taneous cooling and heating ([Reale. Bocchino fc Pered 
l2002t[L6pez-Santiago et al.ll2010[) . The sustained heating 
produces more hot plasma, which then cools and thereby 
enhance the cooler region of the emission measure distri- 
bution. The short, hot flares probably occur in single or 
few loops, and being short, do not ultimately provide as 
large a volume of plasma in the cooler state. 

We should note again the limitations of the three- 
temperature fit: the short flare spectra do have emission 
from lower temperature ions such as O vii, O Viii, and 
Fe XVII (see Figure[7]), and so there should be some fourth 
point in the Figure [8] "Short" curve at low temperatures 
(2-3 MK) with significant emission measure. Such will 
be addressed in future work through detailed emission 
measure reconstruction. 

6. FE K FLUORESCENCE 

At present only the solar corona can be studied in spa- 
tial detail. Fluorescent lines provide a powerful spec- 
troscopic method for probing flare geometry for the un- 
resolved coronae of more distant stars. X-rays of energy 

> 7 keV emitted from a hot corona incident on the under- 
lying photosphere are predominantly destroyed by photo- 
absorption events through inner-shell ionization of atoms 
or weakly ionized species. Observable fluorescent lines 
arise from photons emitted in outward directions by the 
2p— Is decay of the excited atom. In a plasma with solar- 
like composi tion, only Fe fluorescent lines are expected t o 
be detected (|Bailll979l : [Drake. Ercolano fc Swartzll2008| ) . 
The fluorescent line equivalent width depends on the 
height of the X-ray source above the photosphere and 
the heliocentric angle between the flare and line-of-sight. 
Line strength also depends on the relative photospheric 
Fe abundance but is independent of global metallicity ex- 
cept for very metal-poor stars. For coronal scale heights 

> 0.5i?* the Fe K line becomes very weak owing to ge- 
ometric dilution and the larger mean angle of incidence. 
The Fe K line is thus potentially a very powerful diagnos- 
tic of flare geometry and location; it can both constrain 
and calibrate physical models of flares. 

In the flrst detailed study of this kind, Fe K fluores- 
cence was produced by a single large flare in the HETGS 
observa tion of the active giant HR 9024. iTesta et al.l 
(|2007bf) found the origin compatible with photospheric 
photospheric heights and inferred a very compact scale 
height of < 0.1 i?^. A prominent Fe K fluorescent line was 
recently observed from a super flare on the RS CVn-like 



TGCat: EV Lac X-ray Flares 



7 




7^ 



Wavelength (A) 
-a 



» 



o 
o 



13 > 
> X 
O £ 



! XXX 



III I II I II I I I I I I I I I Mil 




o o 
0^ 



1 1 1 1 ■ 

; > 


a 
> 


1 ■ 1 1 1 1 1 1 1 ■ 1 1 1 1 1 1 1 1 


IVII . 

ivn . 


1 1 1 1 ■ ■ 

> > ! 


o 

1 


o 




^ u C " 
1 III 


o o 
1 1 




f 









18.5 



19 



19.5 



20 



20.5 



21 



21.5 



22 



Wavelength (A) 



Figure 7. Photon flux density spectra for different states in a few wavelength ranges: short flares (upper curves; blue), long flares (middle 
curve; red), and the low phase (lowest curve; black). These are "unfolded" spectra, in which the effective area has been factored out, but 
the instrumental line-spread-function remains. The short flare curve has more noise since it has the shortest exposure of 7ks. The others 
have 42 ks (long flares) and 50 ks (quiescent) exposures. Positions of some features (not all actually present) have been marked. The flare 
spectra have been integrated over their durations. 

and hence efficient producer s of Fe K fluorescence from 
theunderlying photosphere (jErcolano. Barlow fc StorevI 
|2005[ ). To search for this fluorescence signature, we 
combined spectra from multiple short flares. The 90% 
confidence upper limit to the Fe K flux is about 6 x 
10^^ phot cm^^s^^. Given our estimate of the 7-20 keV 
flux and the continuum intensity in the vicinity of the 
Fe K line, the expected equiv alent width is ~ 150 eV for 
compact flares computed by iDrake. Ercolano fc Swartd 



binary II Peg by lOsten et al.l (|2007D who found the line 
to be surpri singly strong for photospheric fluorescence. 
lOstenI (|201Cl[ ) detected a very large flare on EV Lac with 
Swift, with a peak Fe Ka flux of 0.027 photons cm~^s~^. 
Their analysis obtained consistent results from fluores- 
cence and hydrodynamic models, with a compact loop of 
height about 0.2R^. 

Based on their rapid decay, the short flares seen from 
EV Lac are expected to be small, single-loop events, 
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Figure 8. The three-component APED fit parameters show the 
changes in temperature (x-axis) and emission measure (i;-axis) for 
the short flares (circles), long flares (stars) and low phase (squares) 
spectra. The short flares are clearly much hotter than the long. 

(|2008f ). corresponding to an Fe K flux of about 2 x 
10~^ phot cni~^s~^. Figure [5] shows our fit to the flar- 
ing spectrum. A more sensitive observation is required 
in order to provide useful constraints from fluorescence 
on the flare scale height. We estimate for double the 
exposure in short flares (14 ks), we can detect a flux of 
2 X 10~^ photcm~^s~^ at 90% confidence. 



HEG; short flares 

6.8 ks r~ 


e XX 


Fe K 










1.6 1.8 2 2.2 


2.4 



Wavelength (A) 

Figure 9. Our current composite flare peak spectrum for 7 ks 
of short flares (dark histogram with light error bars) and best 
fit model (light curve). The best fit to Fe K gives a flux of 
2 X 10~^ photons cm~^ s~^ , and a 90% confldcnce upper limit of 
6 X 10~^ photons cm~^ s~^. 



7. CONCLUSIONS 

We have shown that EV Lac is a reliable factory of X- 
ray flares, and that they come in a large range of shapes 
and sizes, from time scales of ^ 100 s to over 10^ s, and 
a similar range in integrated HETG counts. However, 
their typical brevity and amplitude precludes study of 
spectral evolution for any single flare, so we have re- 
sorted to modeling of composite flares. The evolution 
of temperature and emission measure obtained from line 
emission does indicate that there is sustained heating 
and that simple loop model scaling laws are probably in- 
appropriate. The short flares seem qualitatively different 
from the long flares in their highest weight for the hottest 
plasma. If we assume that there is no sustained heating 
in these flares, then the simple scaling laws would imply 
longer loops due to their higher temperatures, but the 
short timescale imply shorter loops. The simple scaling 
law (see ^ gives Lg ~ 3 — 21, roughly 0.1 to 1 stellar 



radius, a plausible range. Given the diversity of flare 
shapes and scales, we estimate that it would take an ad- 
ditional 200 ks observation to provide sufficient statistics 
in similar moderate length flares for detailed hydrody- 
namic modeling of spectral evolution, and also enough 
counts in short flares to provide a positive detection of 
Fe K fluorescence, and thus provide two independent de- 
terminations of loop sizes. 
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